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ABSTRACT 

Context. Understanding the 3D structure of coronal mass ejections (CMEs) is crucial for understanding the nature and origin of solar 
eruptions. However, owing to the optical thinness of the solar corona we can only observe the line of sight integrated emission. As a 
consequence the resulting projection effects hide the true 3D structure of CMEs. To derive information on the 3D structure of CMEs 
from white-light (total and polarized brightness) images, the polarization ratio technique is widely used. The soon-to-be-launched 
METIS coronagraph on board Solar Orbiter will use this technique to produce new polarimetric images. 

Aims. This work considers the application of the polarization ratio technique to synthetic CME observations from METIS. In particu¬ 
lar we determine the accuracy at which the position of the centre of mass, direction and speed of propagation, and the column density 
of the CME can be determined along the line of sight. 

Methods. We perform a 3D MHD simulation of a flux rope ejection where a CME is produced. From the simulation we (i) synthesize 
the corresponding METIS white-light (total and polarized brightness) images and (ii) apply the polarization ratio technique to these 
synthesized images and compare the results with the known density distribution from the MHD simulation. In addition, we use recent 
results that consider how the position of a single blob of plasma is measured depending on its projected position in the plane of the 
sky. From this we can interpret the results of the polarization ratio technique and give an estimation of the error associated with 
derived parameters. 

Results. We find that the polarization ratio technique reproduces with high accuracy the position of the centre of mass along the 
line of sight. However, some errors are inherently associated with this determination. The polarization ratio technique also allows 
information to be derived on the real 3D direction of propagation of the CME. The determination of this is of fundamental importance 
for future space weather forecasting. In addition, we find that the column density derived from white-light images is accurate and 
we propose an improved technique where the combined use of the polarization ratio technique and white-light images minimizes the 
error in the estimation of column densities. Moreover, by applying the comparison to a set of snapshots of the simulation we can also 
assess the errors related to the trajectory and the expansion of the CME. 

Conclusions. Our method allows us to thoroughly test the performance of the polarization ratio technique and allows a determination 
of the errors associated with it, which means that it can be used to quantify the results from the analysis of the forthcoming METIS 
observations in white light (total and polarized brightness). Finally, we describe a satellite observing configuration relative to the 
Earth that can allow the technique to be efficiently used for space weather predictions. 

Key words. Sun: coronal mass ejections (CMEs) - Sun: corona - techniques: polarimetric - magnetohydrodynamics (MHD) - Sun: 
filaments, prominences 


1. Introduction 


Coronal mass ejections (CMEs) are violent ejections of plasma 
and magnetic flux from the solar corona. They exhibit a wide di¬ 
versity of shapes and properties and there are many unanswered 
questions on their origin and impact on space weather. White- 
light observations from coronagraphs are the most common way 
to detect and study CMEs. However, their 3D structure is not 
fully understood because the solar corona is optically thin. As 
a consequence we observe the 2D projection of CMEs onto the 
plane of the sky. Many previous studies have tried to overcome 
this difficulty by inferring either the density distribution along 
the line of sight (LOS) or the location from where the CME 
emission originates. [Thernisien et al.| ( |2011| ) present a com¬ 
plete survey of the most commonly applied techniques, which 
include forward modelling ( [Thernisien et al.||2009|), polar imet- 
ric ( [Moran & Davila 2()04|), sp ectroscopic ( |Raymond|2QQ2| ), and 
direct inversion ( [Erazin et ar][2009[ ) techniques. More recently, 
unique information on the 3D structure of CMEs have been de¬ 


rived thanks to the observations acquired by the twin STEREO 
spacecraft. In this paper we focus on the v alidation of the polari¬ 
metric technique ( [Moran & Davila|2004[ ). This technique takes 
advantage of the fact that the total and polarized brightness in 
white light from Thomson scattering has a different dependency 
on the angle between the incoming scattered beam and the ob¬ 
server. [Moran et"ar]p010| ) have already compared the results of 
the polarimetric technique with the triangulation t echnique us- 
ing STEREO spacecraft observations. In addition, [Mierla et 
( [201 1[ ) performed an analysis of the limitations of the polariza¬ 
tion technique from an observational point of view. More re¬ 
cently, [Dai et al.] ( [2014[ ) analysed in detail the ambiguities that 
arise in the application of this technique and suggested a method 
for correctly reproducing the CME morphology. 

The coronagraphic images and the stereoscopic observations 
taken from space have proven to be very useful in the reconstruc¬ 
tion of the 3D propagation directions of CMEs. The determina¬ 
tion of this information, as well as the predicted arrival time at 
1 AU, is very important. It is required in order to forecast the 
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possible impact of interplanetary CMEs with the Earth’s mag¬ 
netosphere, thus the possible occurrence of geomagnetic storms. 
Before the availability of STEREO data, many different meth¬ 
ods for inferring the CME propagation direction were developed. 
Examples are the cone model technique applied to LASCO data 
( [Zhao et ST]|2002[ [Michalek et al. [2003 ) and the polarization 
ratio technique ( [Moran & Davila] |2004| ). Since the launch of 
STEREO many other methods have been proposed such as for¬ 
ward modelling ( [Thernisien et al.|2006| ), analysis of height-time 
diagra ms (Mierla et al. 200 8 ), and more recently th e ellipsoid 
model ( [Schreiner et al.|20f3] ). jZuccarello et al.| ( [2^012| ) have used 
the triangulation technique to study the deflection of a CME dur¬ 
ing its path in the solar corona. 

In 2018 the new satellite Solar Orbiter, with the Multi El¬ 
ement Telescope for Imaging and Spectroscopy (METIS) coro- 
nagraph ( jAntonucci et al.|[j6l2 [Eineschi et al.||201^ will be 
able to provide a new view of CMEs. The Solar Orbiter mis¬ 
sion makes it is very important to understand the validity and 
the possible difficulties of the polarization ratio technique for 
the derivation of 3D CME structure and propagation direction. 
Specifically, METIS will observe the coronal emission in polar¬ 
ized white light (bandpass between 580 - 640 nm) and in the UV 
(H I Lyman-alpha 1216 A) at the same time. As Solar Orbiter 
orbit will not be coplanar with the Earth orbit, it will observe 
CMEs from a variety of angles and from different heliocentric 
distances, approaching the Sun to 0.29 AU. As a result the po¬ 
larization ratio technique will be widely used to reconstruct the 
3D morphology of observed CMEs. 

The goal of this work is to apply the polarization ratio tech¬ 
nique to a well-known distribution of density produced from a 
MHD simulation of the ejection of a magnetic flux rope. The 
ejection of magnetic flux ropes is believed to be one of the main 
progenitors of CMEs and this scenario has been confirmed by 
jPagano et al.|p013'^ and jPagano et al.| ( |2013b| ). We first syn¬ 
thesize the white-light total and polarized brightness images of 
the resulting CME as would be seen in the METIS field of view 
when the spacecraft is at its closest approach to the Sun. Next 
we apply the polarization ratio technique to these images and we 
compare the results of the polarization ratio technique with the 
original plasma distribution given by the MHD simulation. This 
will assess the accuracy and the limitations of the technique in 
the 3D reconstruction and trajectory determination of CMEs. Ei- 
nally, we repeat the same procedure for a set of snapshots of the 
MHD simulation to assess the accuracy of the technique in de¬ 
termining the CME properties describing its time evolution, such 
as trajectory, expansion, and velocity. 


2. MHD Simulation 

To simulate the occurrence of a CME and produce a plasma 
density distribution that can assess the validity of the polariza¬ 
tion ratio technique we perform a 3D MHD simulation of a flux 
rope ejection. The plasma density distribution from the resulting 
ejection is used to synthesize METIS observations in white light. 
The simulation is similar to the ones already presented in jPagano 
jet al.| ( [2Q^13a| ), |Pagano et al.| ( |2013b| ), and jPagano et al.| ( |2014| ). 
In these studies a flux rope is first formed through the simulation 
technique of jMackay & van Ballegooij^ ( |2006| ), which is used 
as an initial condition of the MHD simulation. In particular, we 
use the magnetic configuration of Day 19 of the simulation of 
[Mackay & van Ballegooijen| ( |2006| ) when a flux rope has formed 


and is about to erupt. Section 2.2 of jPagano et al. j ( |2() 13 aj ) de 
scribes in detail how the magnetic field distribution is imported 
from the Global Non-Linear Eorce-Eree Eield (GNLEEE)model 
of Mackay & van Ba llegooij^ ( |2006| ) to the MHD simulations. 
jPagano erai.|^2013b| ) identified the parameter space where the 
flux rope is ejected. Einally, jPagano et ar] ( j2014j ) shows how a re¬ 
alistic density and temperature distribution can be constructed to 
account for the presence of a magnetic flux rope. In the present 
work, the simulation parameters are chosen to generate a CME 
in the MHD simulation. 

We use the MPI-AMRVAC software ( jPorth et ^j2Q14j ), to 
solve the MHD equations where external gravity is included as 
a source term. 


dp 
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where t is time, p is density, v velocity, p thermal pressure, and 
B the magnetic field. The total energy density e is given by 


This paper is the second study addressing this topic: In jBem-j 
jporad & Paganoj ( j2015j ) we describe how the position of a single 
blob of plasma is measured using the polarization technique and 
estimate the error associated with the measurement. The present 
paper is a continuation of that study where we apply the same 
processes to a more realistic scenario using MHD simulations of 
a flux rope ejection. The aim is to show how the polarization 
ratio technique and the features of jBemporad & Paganoj ( j20l3] ) 
can be used to interpret polarization ratio images of CMEs. 

The structure of the paper is as follows: In Sect|^we describe 
the MHD simulation; in Sectj^we illustrate how METIS would 
observe the CME produced in the simulation in white light; in 
Sect|^we analyse the results of the polarization ratio technique 
applied to the MHD simulations; in Sect|^we discuss our results 
and draw conclusions. 


P 1 2 


( 6 ) 


where 7 = 5/3 denotes the ratio of specific heats. The expression 
for solar gravitational acceleration is 


g 


GMo 




(7) 


where G is the gravitational constant, Mq denotes the mass of 
the Sun, r is the radial distance from the centre of the Sun, and 
r is the corresponding unit vector. In order to gain accuracy in 
the description of the thermal pressure, we make use of the mag¬ 
netic field splitting te chnique (jPowell et al.|1999] ), as explained 
in detail in Sec.2.3 of jPagano et al.j ( j2013aj ). 

The boundary conditions are treated with a system of ghost 
cells and match those used in [Mackay & van BallegooijelT 
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( |2006| ). Open boundary conditions are imposed at the outer 
boundary, reflective boundary conditions are set at the 0 bound¬ 
aries, and the 0 boundaries are periodic. The 6 boundary condi¬ 
tion does not allow any plasma or magnetic flux through, while 
the 0 boundary condition allows the plasma and magnetic fleld 
to freely evolve across the boundaries. In our simulations the ex¬ 
panding and propagating flux rope only interacts with the 6 and 0 
boundaries near the end of the simulation, thus they do not affect 
our main results regarding the initiation and propagation of the 
CME. At the lower boundary we impose a fixed boundary con¬ 
dition taken from the first four O-cf) planes of cells derived from 
the GNLFFF model. The computational domain is composed of 
256 X 128 X 128 cells, distributed on a uniform grid. The sim¬ 
ulation domain extends over 3 Rq in the radial direction starting 
from r = Rq. The colatitude, 6, spans from 6 = 30° to ^ = 100° 
and the longitude, 0, spans 90°. The simulated domain extends 
to a larger radial distance than the domain used in |Mackay & van| 
|B allegooij en| ( |2006| ) from which we import the magnetic config¬ 
uration. To define the magnetic fleld for r > 2.5 Rq, we assume 
it to be purely radial (Bg = = 0) where the magnetic flux is 

conserved: 

2.5^ 

Br(r > 2.5 Rq,6,(P) = 5,(2.5/?o, 6>, 0) —. (8) 

We set the maximum value of the magnetic fleld intensity to be 
^max = 63 G in the present simulation. 



Fig. 1. Magnetic field configuration used as the initial condition in all 
the MHD simulations. Red lines represent the flux rope, blue lines the 
arcades, green lines the external magnetic field. The lower boundary 
is coloured according to the polarity of the magnetic field from blue 
(negative) to red (positive) in arbitrary units. 


Figure shows a 3D plot of the initial magnetic configu¬ 
ration, which is imported from Day 19 of the simulation by 
Mackay & van Ballegooij en] ( |2006| ). The flux rope (red lines) 
lies in the 6 direction. The flux rope is close to the point where an 
eruption will occur, as it can no longer be held down by the over- 
lying arcades. The arcades are shown by the blue lines, above 
which lies the external magnetic fleld lines (green lines). Some 
of the external magnetic fleld lines belong to the external arcade 
while others are open. For the plasma properties we prescribe 
a non-isothermal solar corona, with a simplified version of the 
treatment described in |Pagano et al.| ( [MT4] ). To produce the tem¬ 
perature distribution we use the function r(B), 


m) = 


2 + (Be/\B\) 


-2 


{Tout — Tmin) + 


(9) 


where Tout is the value of r(B) when Bq = 0. The parameter 
Tynin determines the minimum allowed value for r(B) and occurs 
where Bq = \B\. While the choice of this ad hoc analytic formula 
may seem strange, it is justified by the fact that in our set up the 
flux rope lies in the 6 direction with a positive Bq. It is the only 
structure with a strong magnetic shear in the initial condition of 
our simulation. The form applied in Eq [fallows us to produce a 
cool dense region at the location of the flux rope (i.e. high Bg). In 
principle it is possible to generalize this temperature distribution 
by replacing the 6 direction with the direction of the flux rope 
axis. However, for the present simulations this would have little 
effect. 

The distribution of thermal pressure is independently speci¬ 
fied using the solution for hydrostatic equilibrium with a uniform 
temperature set equal to Tout, 


P = 


- kBlTout exp 

lump 


M^Giump \ / MipOiump \ 

IToutkaRe / \ ^To^kBr / ’ 


( 10 ) 


where plb is the density at r = 7?© when |5| = 0, yu = 1.31 is the 
average particle mass in the solar corona, rup is the proton mass, 
ks is Boltzmann constant. Finally, the density is simply given by 
the equation of state applied to Eq|^and Eq[T0| 


p prnp 


( 11 ) 


p profile 



radial distance 

Fig. 2. Profiles of Logio(p) along the radial direction above the centre 
of the LHS bipole at t = 0 min. The dashed lines show the equivalent 
profiles for a coronal atmosphere in hydrostatic equilibrium at T = 2 
MK. 


In our simulation we choose Tout = 2 MK, Ty^tn = 0.1 MK, 
and Plb = 3.5x10“^"^ g/crrP. With these values, we obtain the at¬ 
mospheric profile shown in Fig|^ this figure draws a radial cut 
of density (solid line) from the lower boundary to the external 
boundary passing through the centre of the LHS bipole (where 
the flux rope lies). As the radial distance increases, the coronal 
density profile decreases as a result of gravitational stratification. 
Below 1.3 /?© there is an excess of density compared to the hy¬ 
drostatic profile (dashed line). This is due to the presence of the 
flux rope that is denser and cooler than the surroundings, as pre¬ 
scribed by Eq|^ The jump at 2.5 Rq in the profile of density 
is a result of the transition between the magnetic fleld configu¬ 
ration interpolated from the GNLFFF model (r < 2.5 Rq) and 
the perfectly radial magnetic fleld (r > 2.5 Rq). Here a small 
Bq component of the magnetic fleld leads to a visible effect in 
the density profile. The dynamic effect of this jump is negligi¬ 
ble, as the coronal dynamics is dominated by the ejection of the 
flux rope. In particular, the density at this location is negligible 
in comparison with the density of the ejected flux rope (See the 
Appendix in [Pagano et al.|2013b] ) . 
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2.1. Evolution of the MHD simulation 


t=0,00m 


The MHD simulation shows a very similar evolution to that 
found in our previous work. The evolution is illustrated in Figj^ 
where maps of density are shown in the (r - 0) plane, passing 
through the centre of the bipoles. In each image we also super¬ 
impose magnetic field lines. Initially the fiux rope lies near the 
lower boundary and as soon as the system is allowed to evolve 
it is ejected upwards. The cause of the ejection is explained in 
|Pagano et al.| ( |2013a| ) where we describe how the magnetic stress 
built up during the fiux rope formation is released. Once the 
ejection occurs, a high density structure rises and expands. Ini¬ 
tially, the high density region maintains the shape of the ejected 
flux rope (Fig|^); however, near the end of the simulation it is 
less identifiable with the fiux rope (Fig|^). After approximately 
100 minutes, the high density region reaches the outer bound¬ 
ary at 4 7?©. It should be noted that the fiux rope propagates 
non-radially (yellow dashed line in Figs, [^-c) in the direction 
of the null point, which lies above the arcade system. It fol¬ 
lows the non-radial path as the confining Lorentz force exerted 
by the magnetic arcades is weakest there and so this becomes 
a favourable escape direction. This leads to a defiection of the 
flux rope by a few degrees from a perfectly radial propagation. 
The magnetic field configuration undergoes a major evolution 
and reconfiguration as a result of the ejection. Magnetic fiux is 
expelled outwards and at = 69.60 min (Figj^) a region of com¬ 
pressed magnetic field is visible ahead of the ejection, where the 
front of the ejection compresses both plasma and magnetic fiux. 

The fiux rope ejection leads to an increase in density in the 
outer corona that propagates outwards and can be labelled as a 
CME. It also exhibits the three-part density structure associated 
with CMEs. Eigure|^ shows the density profile along the line of 
propagation of the flux rope at each of the times shown in Eig|^ 
The dashed vertical lines show the position of the centre of the 
fiux rope at each time. In particular dX t = 69.60 min we have 
a CME front where the density has increased by approximately 
two orders of magnitude compared to the value of density at the 
same location at ^ = 0. Behind the front there is a void, where 
the density is about half that found in the front. Eurther behind 
the fiux rope constitutes the core of the CME, where the density 
is approximately one to two orders of magnitude higher than that 
of the void. Einally, Eigj^ shows the position of the centre and 
top of the fiux rope as a function of time along the line of prop¬ 
agation. The centre of the fiux rope is given by the maximum 
of the quantity Bel\B\ along the propagation line and the top is 
located where the quantity Bol\B\ changes sign above the centre. 
The parameters we have chosen for the simulation give condi¬ 
tions in which the fiux rope is expelled from the solar corona at 
a typical velocity. The fiux rope undergoes an acceleration phase 
for the first 30 minutes and then it travels outwards at a constant 
speed. The overall average speed of the fiux rope is about -157 
kmjs, whereas its front travels slightly faster at - 274 kmjs, as 
it expands whilst propagating outwards. 


3. Synthesis of METIS observations in white iight 

The synthesis of total white-light (tB) and polarized (pB) bright¬ 
ness images is carried out in 3 steps. In step one we interpo¬ 
late the outputs of the spherical 3D MHD simulation onto a 3D 
cartesian grid. Eor the second step, we compute the total and 
polarized brightness of each plasma element in the observer’s 
direction. Einally, for the third step we integrate along the line 
of sight (LOS). 
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Fig. 3. Maps of Logioip) in the (r - 0) plane passing through the centre 
of the bipoles - 34.80, and t - 69.40min. Superimposed are 

magnetic field lines plotted from the same starting points (green lines). 
The yellow dashed line shows the flux rope propagation line. 
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Fig. 4. Profiles of Logio(p) along the yellow dashed line in Fig[3] 
which represents the direction of the flux rope ejection at ? = 0 min 
(black line), t = 34.80.60 min (blue line), and t = 69.60 min (red line). 
The dashed lines mark the position of the centre of the flux rope at each 
time. 


The first step is carried out through a 3D interpolation that 
has already been successfully applied in [Pagano et al.| ( |2014| ). 
This interpolation is carried out to allow the integration along 
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Fig. 5. Position of the radial distance of the centre and top of the flux 
rope as a function of time. 


the LOS of the originally spherical domain, thus transforming 
the domain of the spherical MHD simulation into a cartesian 
grid. In this grid the flux rope lies on the solar disk at about 
24° from the plane of sky (POS) in the direction of the observer. 
This can be done as the longitudinal coordinate is arbitrary. The 
choice of 24° is made to study a general case where the CME 
is not ejected too close to nor too far away from the POS. The 
resulting cartesian grid has an origin at the centre of the Sun, z is 
the direction along the LOS, y is the direction parallel to the NS 
axis of the Sun, and the x direction is perpendicular to both z and 
y. Owing to geometrical effects, the spherical simulation cannot 
fill every point of the cartesian box with values. Therefore, to 
complete the cartesian box the initial density and temperature 
distribution are mirrored about the planes: 0 = 0°, 0 = 90°, 
6 = 30°, and 6 = 100°. Multiple mirrored copies of the plasma 
distribution for the MHD domain at r = 0 are used to populate 
part of the cartesian box that cannot be sampled from the spheri¬ 
cal simulation. In such a way i) plasma is present in the cartesian 
box everywhere (r < 4 Rq), and ii) outside the domain specified 
by the spherical simulation, the corona remains unperturbed as 
given by the initial condition. 

Figurej^shows a generic LOS (red dashed line) drawn on the 
corresponding density map on the v - z plane at y = 0.3 /?© and 
t = 69.60 min. The flux rope ejection produces many distinct 
density structures travelling outwards and the LOS crosses some 
of them. In contrast, the corona outside of the spherical MHD 
simulation domain is unperturbed and remains in its initial state. 

The second step is to compute the tB and pB emissions from 
each plasma element of the cartesian box. This is carried out by 
computing the Thomson scattered light from a single electron in 
each cell of the cartesian box, depending on the impact distance 
of the LOS and the angle from the plane of sky ( |Minnaert|1930| ). 
The contribution is then multiplied by the electron number den¬ 
sity of each cell. We point out here that the white-light emission 
of a plasma element only depends on the plasma density, its po¬ 
sition with respect to the solar disk, and the specific scattering 
angle. Figure [7 shows the computed tB and pB emissions along 
the LOS in Fig ^ The vast majority of the emission is from the 
region within I Rq from the plane of sky (z = 0). Therefore, 
the assumptions made to populate the peripheral regions of the 
cartesian box have no significant impact on our study. 

The final step consists of integrating the contribution from 
each plasma element along a given LOS. In the present work 
we use 320 x 400 lines of sight and we choose to show only 
the lines of sight crossing the METIS field of view when Solar 
Orbiter reaches perihelion at 0.29 AU from the Sun (I.6/?o < 
< 3.1/?o). 

Figure shows the synthesized observation of METIS in 
white-light from our MHD simulation at t = 69.60 min, where 
the red contour highlights the approximate location of the centre 
of the fiux rope axis, and the yellow contour is the CME front. 
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Fig. 6. Map of Logioip) in the cartesian box at y = 0.3 Rq with a LOS 
diix - 2 Rq overplotted with red dashed line. 
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Fig. 7. Emission in white light tB (black line) and pB (blue line) along 
the LOS shown in Eig[^ 


The two blue lines mark the METIS field of view. At this ra¬ 
dial distance from the Sun, the fiux rope is still well defined and 
visible in the observations as the brightest region of the ejection. 
METIS will observe close enough to the solar surface to catch 
coherent structures before they merge into the solar wind stream 
of the outer corona. In contrast, the front of the CME seems 
quite faint in our case study. Figureshows the polarized com¬ 
ponent of the white-light and this mostly reproduces the same 
structures seen in Figj^, except with lower brightness. Both 
panels in Figj^show a sharp discontinuity in emission at the 6- 
borders of the simulation domain. This is due to the limited 6 
and 0 extension of our spherical domain, as in these regions we 
have assumed the corona to be maintained unperturbed. 
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Fig. 8. (a) Synthetic image of the total white-light brightness {tB) 

integrated along the line of sight in the METIS FOV. The red contour 
roughly identifies the fiux rope position where the ratio between the 
y component of the magnetic field integrated along the line of sight 
divided by the magnetic field intensity is above the arbitrary positive 
threshold of 0.4, and the yellow contour the position of the ejection 
front where this quantity is lower than the arbitrary value of -0.43. (b) 
Synthetic image of the polarized white-light brightness {pB) integrated 
along the line of sight in the METIS FOV. Both images are synthesized 
from the MHD simulation at t = 69.60 min and using a METIS field of 
view at perihelion at 0.29 AU from the Sun and extending from 1.6 Rq 
to 3.1 (blue lines). 


4. Results 

From the synthesis of METIS observations, we obtain tB and pB 
images that are equivalent to actual white-light observations of 
the solar corona with the METIS coronagraph. In this section, 
we present the results obtained from the analysis of synthesized 
images with polarization ratio and excess brightness techniques. 
To derive the CME plasma 3D distribution and electron density, 
we compare the results with the actual density distribution from 
the MHD simulation. Finally, we interpret the comparison in 
terms of the analysis carried out in |Bemporad & Pagano| ( [2Q15] ). 


4.1. Reconstruction of the 3D CME structure and 
propagation direction 


We apply the polarization ratio technique to synthesized METIS 
difference images, i.e. images from which the pre-CME quiet 
corona emission has been subtracted. In observations, this sub¬ 
traction is important in order to isolate the so-called excess 
brightness (i.e. the total brightness in white light due to the 
CME plasma alone tBcM^)- h removes the spurious emission 
due to background coronal plasma located along the same LOS 
but external to the CME volume. In this study, the white-light 
image corresponding to the MHD simulation at r = 0 is used 
as the quiet corona, thus our assumption of the density distribu¬ 
tion outside of the spherical simulation becomes irrelevant for 
the results of the analysis. 

We identify a sample of lines of sight to apply the polariza¬ 
tion ratio technique. Figurej^ shows the identified lines of sight, 
all intercepting our MHD simulation box and showing a white- 
light brightness above a given threshold. These lines of sight 
identify a distinct set of (x,y) points. For each (x,y) point the 
polarization ratio technique gives a corresponding position z of 
the emitting plasma along the LOS, thus a cloud of points (x,y,z). 
This cloud of points is compared with the cloud of points corre¬ 
sponding to the centre of mass of the density distribution folded 
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about the plane of sky (call ed hereafter folded centre of mass) 
( [Bemporad & Pagano|2Q15| ): 


CMfoidedi^^y) — 

z[(p(z, X, y, t) - piz, X, y, 0)) + (p(-z, x, y, t) - p{-z, x, y, 0))] dz 
J^(p(z, X, y, t) - p(z, X, y, 0)) + (p(-z, x, y, t) - p(-z, x, y, 0)) dz 


( 12 ) 


Thus we have a cloud of points (x, y, z) resulting from the appli¬ 
cation of the polarization ratio technique (red crosses hereafter) 
and a second cloud of points (x, y, z) (green crosses hereafter) 
computed from the positions of the folded centre of mass along 
each LOS in the cartesian simulation box. The two clouds are 
constructed in a way to have the same (x, y) values, but can dif¬ 
fer in z values. Figure shows the position of these points as 
seen by the observer. In Fig|^ we show the two clouds of points 
projected onto the x - z plane (corresponding to a view from 
the north pole), and in Fig|^ we show the two clouds of points 
projected onto the z - y plane (corresponding to a view perpen¬ 
dicular to the equatorial axis and the observers point of view). 
In ah plots, red crosses represent the cloud of points obtained 
from the polarization ratio and the green crosses represent the 
cloud of points obtai ned from Eq[^ As explained in |Bempo-| 
rad & Pagano| ( |2Q15) ), the position of the folded centre of mass 
is the best approximation of the position output obtained from 
the polarization technique. This is mainly due to the ambigu¬ 
ity in the polarization ratio expression where it is not possible 
to distinguish whether the emitting plasma is located in front or 
behind the plane of the sky (POS). In the following discussion 
we estimate the error associated with this approximation. 

The polarization ratio technique catches the 3D cloud of 
points of the folded centre of mass with sufficient accuracy to 
describe the macroscopic structure, as the two clouds of points 
mostly overlap. Also, the shapes are nearly identical from both 
points of view perpendicular to the LOS (Fig|^, Fig|^). More¬ 
over, from Fig|^ we see that both clouds of points are aligned 
along a direction of 25° from the plane of sky below r = 1.6 Rq 
and the clouds become more dispersed around a direction of 15° 
from the plane of sky. As can be seen in Fig|^ and Fig|^, the 
two clouds are nearly identical in shape and position. However, 
there is an offset between the two clouds, as the cloud of points 
obtained from the polarization ratio (red crosses) is about 2° far¬ 
ther from the POS than the cloud obtained from the folded centre 
of mass. Figureshows that the two clouds follow this pattern 
independently of latitude. 

A more detailed comparison can be carried out to measure 
the errors. We consider a smaller sample of lines of sight, ah 
lying on the same y coordinate. Figure shows a set of points 
taken at y = 0.9Rq (crossing the centre of the flux rope) and 
in the image we indicate the respective angle with the plane of 
sky for a sample of representative points. Green and red crosses 
never depart from one another by more than O.l Rq where the 
distance is the greatest at the outer boundary of the METIS held 
of view. At the outer boundary the angle with the plane of sky 
is lower than 20°, as explained in [Bemporad & Pagano| ( |2015) . 
This leads to an overestimation of the distance from the plane 
of sky of the folded centre of mass position (red crosses are far¬ 
ther from the plane of sky than green crosses). Instead, at the 
location of the flux rope (brightest features) the agreement be¬ 
tween the points obtained by the polarization technique and the 
folded centre of mass is much better. In this region the struc¬ 
tures he at an angle of ~ 20° or more with the plane of the sky. 
As described in [Bemporad & Pagano| ( |2015] ) the polarization ra¬ 
tio technique slightly underestimates the distance from the plane 
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Fig. 9. Maps of Log\Q of column density (range between 10"^^ g/cm^ 
and 10“^^ glcm^) seen from three different points of view; superim¬ 
posed is the cloud of points inferred from the polarization ratio tech¬ 
nique (red crosses) and the cloud of point computed from Eq[T^ (green 
crosses) projected onto the plane of sky: (a) point of view parallel to 
direction z and (v,y) coordinates of clouds of points; (b) point of view 
parallel to direction y and (v, z) coordinates of clouds of points; (c) point 
of view parallel to direction v and (z,y) coordinates of clouds of points, 
(d) Maps of white-light emission seen from a point of view parallel to 
the y direction; superimposed is a subset of points from the two clouds 
at y = 0.9 Rq. In light blue we give the angle with the plane of sky for 
some of the points. The thick magenta line marks the plane of sky. 


of sky in such circumstances (green crosses are farther from the 
plane of sky than red crosses). However, this effect is not as evi¬ 
dent as in |Bemporad & Pagano| ( |20I5] ) owing to the presence of 
multiple structures along a line of sight, which marks a signifi¬ 
cant difference from the idealized situation of a single blob. 


Figure [T^ shows the differences in the z coordinate between 
the two clouds of points. The differences in the position can be 
up to ~ 0.1 7?©; however, at the location of the flux rope we find 
that the differences are much smaller ( < 0.03 Rq). This is due to 
the increased accuracy of the polarization ratio technique when 
the source of the white-light emission is highly localized. Over 
the whole field of view, the polarization ratio technique tends to 
overestimate the calculated distance from the plane of sky com¬ 
pared to the position of the folded centre of mass. In contrast it 
should be noted that the distance from the plane of sky is under¬ 
estimated in the region where the flux rope is present. This oc¬ 
curs because in this region the folded centre of mass lies at more 
than 20° from the POS (see Fig|^ and Fig|^) and the measure¬ 
ment from the polarization ratio leads to an underestimation of 
the distance from the POS. 


4.2. Comparison with the centre of mass 

In Sect |4.1| we illustrate that the polarization ratio technique re¬ 
constructs the position of the folded centre of mass, but intro¬ 
duces a measurement error of the order of ~ 0.03 Rq. How- 
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Fig. 10. Map of differences of z coordinates in Rq for the two clouds 
of points, overplotted on a map of Logio of column density as in Fig[^. 
Red regions correspond to where the distance from the plane of sky 
computed from the polarization ratio technique is overestimated and 
blue regions correspond to where it is underestimated. The magenta 
contour roughly identifies the flux rope position where the ratio between 
the y component of the magnetic field integrated along the line of sight 
divided by the magnetic field intensity is above the arbitrary positive 
threshold of 0.4. 


ever, an outstanding question is how these measurements com¬ 
pare with the true position of the CME, i.e. the centre of mass 
along the LOS. 

In Fig[^we show the images equivalent to Fig|^-c where 
we compare the cloud of points obtained by the polarization ratio 
technique (red crosses) with the cloud of points determined from 
the centre of mass (light blue points). The two clouds of points 
in Figshow a very similar shape and extension; however, they 
have a systematic offset, significantly larger than the one high¬ 
lighted in Fig|^ In particular, the polarization ratio technique 
systematically overestimates the distance from the plane of sky 
of the centre of mass and the overestimation seems greater at 
larger radial distance from the Sun (FigfTT^). For instance, this 
is evident for the cloud points below 1.5 7?© that correspond to 
the location of the flux rope. The light blue points (centre of 
mass) lie between 20° and 25°, which is where we have posi¬ 
tioned the flux rope in our cartesian box. In contrast, the red 
points (output of the polarization technique) lie between 25° and 
30°, hence the distance from the POS is overestimated. 

No clear dependence on latitude is seen (Fig[TTJ)). In this 
study the two clouds depart by about 10°, but this value strongly 
depends on the amount of CME plasma that is present beyond 
the POS. This is a direct consequence of the initial ejection po¬ 
sition of the flux rope as well as its expansion. In general the 
matching between the two clouds is significantly worse than in 
Fig 1^ and the offset between the clouds in Fig[TT]is ~ 0.2 Rq. 


4.3. Electron density estimate and comparison with MHD 
simulation 

Another application of white-light images is to measure the den¬ 
sity of the plasma of the CME. The electron density has been 
computed (as is usually done for real observations) from the ex¬ 
cess brightness tBcM^ observed inside the simulated CME. This 
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Fig. 11. Maps of Logio of column density (range between 10"^^ g/cm^ 
and 10"^^ glcm^) seen from two different points of view; superimposed 
are the cloud of points inferred from the polarization ratio (red crosses) 
and the cloud of points of the centre of mass along the LOS (light blue 
crosses) projected on the plane of sky: (a) point of view parallel to the y 
direction and (v, z) coordinates of the clouds of points; (b) point of view 
parallel to the x direction and (z,y) coordinates of the clouds of points. 


quantity is given by an integral of the white-light brightness due 
to ah the electrons aligned along the considered LOS; the elec¬ 
tron column density Ne (cm“^), corresponding to the LOS inte¬ 
gral of the electron density distribution Ue (cm“^), is thus pro¬ 
vided by dividing the IBcme by the brightness of a single elec¬ 
tron tBg. However, owing to the geometry of Thomson scatter¬ 
ing, the latter quantity depends on the location of the considered 
electron along the LOS; in particular, tB^ maximizes at the POS 
and then progressively decreases when moving away from the 
POS along the LOS. The usual assumption in the data analysis 
is that the bulk of the emission comes from the plasma located 
around the POS (as also shown in Fig|7]). Therefore a lower 
limit estimate of the column density Ne is obtained by dividing 
the excess brightness by the tB^ exactly on the POS (z = 0). 
We now explore a possible rehnement of this technique: we es¬ 
timate the column density Ne by dividing the excess brightness 
tBCM£ by the single electron brightness tB^ assuming that it is 
located on the position along the LOS determined from the po¬ 
larization ratio technique. Hence, we have two slightly different 
measurements of the column density. First, by assuming that 
the white-light emission originates from a single bulk of plasma 
placed either in the POS or second, at the position along the LOS 
computed by the polarization ratio technique. Such assumptions 
are reasonable in order to achieve an estimation of the plasma 
involved in a CME. However, we need to verify that they can 
reproduce an actual map of the column density in a given held 
of view. 
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Fig. 12. (a) Maps of Logio of column density number (A^^) computed 

from the MHD simulation, (b) Maps of Logio of column density num¬ 
ber (Ne) computed from the technique applied to the MHD simulation. 
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In FigfT^ we show the actual column density computed 
from the MHD simulation and in Fig[T^ we show the column 
density computed from synthesized METIS images assuming the 
plasma is ah placed in the plane of sky. The two column den¬ 
sity distributions show a high degree of similarity over the range 
of densities investigated here. In general the magnitude of the 
column density is reproduced, as well as its distribution and fea¬ 
tures. 



Fig. 13. Maps of the contrast (Eqp^ between column density 
computed by synthesis of white-light images and the actual column den¬ 
sity from the MHD simulation under two different assumptions: (a) all 
emitting plasma is located on the plane of sky; (b) all emitting plasma is 
located at the distance from the plane of sky computed by the polariza¬ 
tion ratio technique. The white contour roughly identifies the flux rope 
position where the ratio between the y component of the magnetic held 
integrated along the line of sight divided by the magnetic held intensity 
is above the arbitrary positive threshold of 0.4. 


In order to illustrate the error associated with this technique, 
we map the contrast between the column density computed by 
the polarization ratio technique and the actual column density 
computed from the MHD simulation in Fig[^ The contrast 
is computed as 


Ne(WL) - Ne(MHD) 
Ne(MHD) 


(13) 


where Ne(WL) is the column density computed from the synthe¬ 
sis of white-light images and Ne(MHD) is the column density 
of the MHD simulation. We consider the quantity Ne(WL) un¬ 
der two assumptions: i) the plasma is centred on the plane of 
sky (Fig[T^, “POS assumption”) and ii) the plasma is centred 
on the position along the LOS as computed by the polarization 
ratio technique (FigfT^, “LOS assumption”). In both maps we 
also show the position of the flux rope (white contour). 

The distribution of the error shows some interesting features. 
We find that the POS assumption leads to an error of up to 10% 
on the column density value. In the region above the flux rope 
position the column density is generally underestimated by both 
techniques by about 5%. At the location of the flux rope the 
underestimation is around 10% for the POS assumption, but re¬ 
duces to ~ 3% in the LOS assumption. The column density is 
only overestimated in a small location near the equator under 
the LOS assumption. The underestimation near the flux rope is 
larger for the POS assumption (Fig[T^) than for the LOS as¬ 
sumption (FigfT^) because the approximation of the latter (a 
bulk of plasma at the position determined by the polarization ra¬ 
tio) is more precise than the approximation of the former (a bulk 
of plasma at the plane of sky). The same is true for all the loca¬ 
tions that presented a dense feature in Fig[T^ as in general the 
relative error associated with the LOS assumption is half of 
that associated with the POS assumption at any given location. 
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This interesting result suggests an improved technique for the 
determination of column densities from white-light images of 
CMEs. In particular, uncertainties can be minimized by assum¬ 
ing that in each pixel of the 2D coronagraphic image the emitting 
plasma is not centred on the POS (as is usually assumed in the 
literature), but is centred on the actual location along the LOS as 
inferred from the polarization ratio technique. 

4.4. Time evolution 

By taking into account a series of snapshots from the MHD sim¬ 
ulation, processing t hem a s in Se c|3]an d applying the same anal¬ 
ysis described in Sec |4.I| and Sec |4.3| we can investigate the ca¬ 
pabilities of the polarization technique in inferring the evolution 
of the CME. In particular we can consider its trajectory, velocity, 
expansion, and the mass involved as a function of time. 

In contrast to Sec|^ we now want to apply the analysis to 
lines of sight that cross the CME, and to do so we need a way 
to identify them. We decided this was the best choice because 
the polarization ratio technique performs more accurately when 
applied to lines of sight that cross dense structures, which are 
the most relevant for our study. The polarization ratio pB/tB of 
the white-light emission from a single structure is highest on the 
POS and decreases with distance from the POS. In our simula¬ 
tion the CME is composed of structures that do not he on the 
POS. We use this property to identify the CME in the LOS inte¬ 
grated images (panels in Eig. synthesized from the snapshots 
of the MHD simulation. The CME cloud is identihed as the re¬ 
gion where the polarization ratio of the LOS integrated images is 
less than a critical value of pB/tB = 0.6. It can be shown that be¬ 
yond a certain heliospheric distance, the ratio pB/tB of the emis- 
sion decreases while a structure moves radially (Eig.6 in |Bem-| 
Iporad & Pagano| ( [20I5| )). Therefore the use of a constant critical 
upper value for pB/tB is a reasonable technique for identifying 
structures that are mostly radially propagating, like CMEs. A 
map of polarization ratio with the considered lines of sight at 
t = 127.6 min is shown in Eig[^ Erom this point of view, the 



X [Rq] 

Fig. 14. Map of polarization ratio (pB/tB); superimposed are the lines 
of sight that interject the CME according to our criteria {pBjtB < 0.6) 
at t = 127.6 min. 

CME appears as an expanding cloud wider near the solar surface 
with a cusp that represents the farthest distance reached by the 
CME. This method is applied ad hoc for this case in order to 


describe a smooth CME evolution and remove spurious features 
from the analysis. It does not aim at being a general method 
for CME detection in coronagraph images. Different approaches 
are possible, and more complex detection techniques exist (e.g. 
Bonte et al.|2011| ). Erom the polarization ratio we know that the 
CME enters the METIS held of view after only ~ 80 min and 
stays mostly inside the held of view for an additional 60 min. 
Consequently we consider in our analysis only snapshots within 
this time range. 

Eigure shows the CME cloud of points at different times 
from three different points of view (POS, north pole, west limb). 
We use a different colour intensity for each cloud to represent 
different snapshots (from dark to bright), so that the time evolu¬ 
tion can be perceived. Eigure p~5^ shows that our CME clouds 
are composed of a compact front when the CME just enters the 
METIS held of view and it evolves into a less compact front with 
a triangular shape. The time evolution of the three clouds gen¬ 
erally conhrm the analysis carried out in Sec|^ Erom the Sun’s 
north pole view (EigfT5|)) we see that the CME keeps a narrow 
shape along the z-direction and it propagates at farther radial dis¬ 
tances with a coherent motion. Eigure \T5\ ) also shows that the 
clouds of points from the folded centre of mass and from the po¬ 
larization ratio technique mostly coincide with a small offset that 
displaces the folded centre of mass slightly farther from the POS 
than the polarization technique. As we found in Pig[T^ the lines 
of sight that cross the flux rope show an inverted behaviour with 
regard to the other lines of sight. Namely that the folded centre 
of mass is farther from the POS than the position obtained from 
the polarization technique. This is the reason why the relative 
position between the red and green clouds is switched between 
Pig [T^ and Pig|^ This happens as the flux rope structure is sig- 
nihcantly denser than the background corona and it lies at an 
angle larger than 20° from the POS than is shown in Pig|^. 

At the same time, both clouds (green and red) are displaced 
signihcantly farther from the POS than the true centre of mass 
(light blue). Pigur^T^ shows again the offset along the z direc¬ 
tion and also that the three clouds do not present any signihcant 
discrepancy along y. It also reveals the rough bow shape of the 
propagation front of the CME. 

Similarly, we plot in Pig[T^the angle from the POS as a func¬ 
tion of the heliocentric distance projected onto the POS in order 
to show the progression of the CME and the accuracy of the po¬ 
larization technique in catching its trajectory. In order to trace a 
rough trajectory of the CME, we represent each snapshot by the 
LOS with the maximum projected radial distance and its angle 
with the POS (red, green, and light blue lines in Pig[T^. We hnd 
that all the clouds follow the same pattern, where the CME is de¬ 
flected by ~ 5°. The trajectory of the cloud of the centre of mass 
(light blue) starts at ~ 24° (where the flux rope is placed) and it 
ends at ~ 29°. The trajectories of the clouds of the polarization 
ratio technique and folded centre of mass show an offset of about 
2° and 3°, respectively, from this position. Prom this estimate of 
the trajectory, we can infer a CME velocity of about 110 km Is 
and a velocity projected onto the POS of 130 km/s. This has to 
be compared with the velocity of the front of the CME as com¬ 
puted in Pig|^ where we found that the speed of the flux rope 
was of 157 km/s. 

Eigure shows the angular extension of the CME for the 
three clouds as a function of time where we also report the ex¬ 
pansion factor at given times. The expansion factor is simply the 
ratio between the extension of the cloud and its initial extension. 
We And that the centre of mass cloud and the folded centre of 
mass cloud always show a similar extension, while the polariza¬ 
tion technique cloud is always narrower. However, this does not 
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Fig. 15. Maps of LogiQ of column density (range between 10"^^ 
glcrn^ and 10"^^ glcm^) seen from three different points of view; super¬ 
imposed are the points of the clouds that belong to the CME. Points ob¬ 
tained from the polarization ratio technique are red, points of the folded 
centre of mass are green, and points of the centre of mass are light blue, 
(a) Point of view parallel to direction z and (x, y) coordinates of clouds 
of points; (b) point of view parallel to direction y and (x, z) coordinates 
of clouds of points; (c) point of view parallel to direction x and (z,y) 
coordinates of clouds of points. A movie of these panels is available 
online. 


significantly affect the measurement of the expansion, which in 
the end differs by only 0.1. The folded centre of mass and centre 
of mass clouds present an expansion about 5% larger than the 
polarization technique cloud. 
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Fig. 16. Angle from the POS as a function of the projected heliospheric 
distance for the clouds of points as in FigffS] 
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Fig. 17. Angular extension of the CME as a function of time accord¬ 
ing to the different clouds. Superimposed are the expansion factors at 
certain times. 
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Fig. 18. Average column density per line of sight as a function of time. 


Finally, we can sum the column density measured by the po¬ 
larization technique over the CME lines of sight and compare it 
with the actual column density along the same lines of sight in 
the MHD simulation. This will allow us to assess the mass vari¬ 
ation of the CME as a function of time. During its propagation, 
as the CME expands in the METIS FOV, it covers a wider area 
and the total column density involved in the CME at the end of 
the simulation is 5 times greater than the mass involved when it 
enters the METIS FOV. Since the mass measurement is greatly 
affected by the number of lines of sight that we consider belong 
to the CME, Fig[^ shows the total column density of the entire 
CME as function of time where it is divided by the number of 
CME points at each time. It shows the value computed from the 
polarization technique under the LOS assumption (green line) 
and from the MHD simulation (red line). The averaged column 
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density steadily decreases and by the end it is ~ 92% less than 
the value found at the beginning. However, the discrepancy be¬ 
tween the two measurements is small where its maximum value 
corresponds to only ~ 3% of the average column density. This 
proves once more that the polarization technique is a reliable 
technique for measuring the column density along a line of sight. 


5. Discussion and conciusions 

In this work we have investigated the reliability of results of 
the 3D structure of CMEs as inferred from the polarization ra¬ 
tio technique. The technique was applied to white-light images 
of total (tB) and polarized (pB) brightness synthesized from a 
MHD simulation of a flux rope eruption. For synthesis purposes 
we used the METIS coronagraph held of view. By comparing the 
real plasma density distribution with the results from the polar¬ 
ization ratio technique, we analyse how the position, direction of 
propagation of the CME, and the column density of plasma are 
reconstructed. 


5.1. Observations of CMEs by METIS/Solar Orbiter 

The synthesis of the white-light images as seen by METIS from 
the MHD simulation allows us to show possible future obser¬ 
vations of CMEs when the Solar Orbiter spacecraft will reach 
perihelion at 0.29 AU. At this heliocentric distance the projected 
held of view of METIS will extend down to 1.8 Rq from the 
solar surface. This will observe an altitude where in general 
ejected flux ropes have not yet merged in the solar wind. Hence, 
METIS will allow us to carry out coronagraphic observations of 
fully formed ejected flux ropes in both visible and UV ranges. 
Furthermore, both METIS (imaging the solar corona) and EUI 
(imaging the solar disk) will observe with the Ly-a intensity. 
Solar Orbiter will therefore provide a unique opportunity to fol¬ 
low the evolution of flux ropes from their formation to ejection 
out to about 3 Rq where the merging of the flux rope with the 
solar wind eventually occurs. At the same time, the METIS held 
of view will be large enough to catch in the same image the 
CME front and other related phenomena such as CME-driven 
shocks and side interactions - reconnection with nearby coronal 
structures. With the spatial resolution of the instrument (20 arc- 
sec/pixel) enhanced by the proximity of the satellite to the Sun, it 
will be possible to describe with unprecedented clarity the den¬ 
sity distribution of the three components of CMEs. In addition 
we will be able to derive information on the CME plasma tem¬ 
perature by combining for example white-light and UV images 
acquired at the same time. In particular, this latter point will be 
the subject of future work where synthetic white-light and Ly-a 
METIS images of CMEs will be analysed in order to enhance 
our diagnostic capacity and to develop techniques that will be 
available for when actual data is available. In the work presented 
here, we use a MHD simulation and the synthesis of visible-light 
images to show preliminary aspects of how METIS will observe 
CMEs when in operation. Additionally, it is worth considering 
that Solar Orbiter will travel on an orbit not co-planar with the 
Earth’s orbit. Although we have carried out our analysis from a 
point of view on the Earth’s orbit, it is possible to imagine that 
the polarization technique can be used to study the 3D structure 
of CMEs from points of view off of Earth’s orbit plane. 


5.2. 3D Reconstruction of a CME with polarization ratio 


Overall, the results from the polarization ratio technique suc¬ 
cessfully infer the 3D CME structure once we have a clear under¬ 
standing of the interpretation of results and the errors associated. 
When applied to the present MHD simulation, the polarization 
ratio technique is able to infer the location of the flux rope and 
other density structures with sufhcient precision to describe the 
general structure of the CME. In particular we confirm here, as 
already discussed in |Bemporad & Pagano| ( |2015] ), that the polar¬ 
ization ratio technique provides a good estimate for the location 
along the line of sight (LOS) of what we call the folded centre of 
mass, namely, the centre of mass of the density distribution ob¬ 
tained by summing the distribution behind the plane of the sky 
(POS) to the one in front of it. In particular, the technique mea¬ 
sures the position of the folded centre of mass within an error 
of ~ 0.03 Rq. This is a relatively small error considering that 
the extension of a CME can be over I Rq in the radial direction. 
Moreover, considering a CME velocity on the order of 500 km/s, 
a distance of ~ 0.03 Rq is covered in approximately 40 seconds, 
while the CME event lasts up to a few hours. Thus, from this 
point of view, the use of the polarization technique would lead 
to a minimal error if applied to the determination of large-scale 
properties such as the injection time of the CME into the solar 
wind or the arrival time at the Earth’s magnetosphere. 

Nevertheless, it should be noted that such a precision is re¬ 
alistic only when the folded centre of mass and the centre of 
mass are placed approximately at the same location. This only 
occurs when the CME travels and expands without crossing the 
POS. If, however, the CME bubble crosses the POS the errors 
introduced by the polarization ratio technique may be up to ~ 10 
times larger. This is shown in Fig [T9] where the CME bubble 
(assumed to be spherical and with homogeneous density) lies 
entirely out of the POS (top left) the location of the real cen¬ 
tre of mass (red dashed line) is coincident with the centre of 
mass of the folded density distribution and with the one inferred 
from the polarization ratio technique (orange dashed line). Nev¬ 
ertheless, when the CME bubble partially crosses the POS (top 
right) a disagreement between the two distributions starts to ap¬ 
pear. The disagreement is larger at the centre of the CME where 
the observed tB and pB intensities come from a longer integra¬ 
tion along the LOS. The uncertainty maximizes when the CME 
expands along the POS (bottom left); due to the symmetry in 
Thomson scattering geometry, it is not possible to distinguish the 
case when the CME is expanding in front of (top right) or behind 
(bottom right) the POS. It should be noted here that the estima¬ 
tion of the errors associated with the polarization ratio technique 
that we And here, using a realistic MHD simulation are in agree¬ 
ment with those derived in [BemporaT^Taga^ ([20^. This 
is an important cross-check that validates our previous analysis. 
|Dai et al.j ( |20I4| ) have also concluded that the polarization ratio 
technique is affected by an ambiguity in identifying the location 
of the source of emission and have proposed a method to over¬ 
come this. Our work does not tackle the problem from this gen¬ 
eral point of view, but we focus on the errors that persist in the 
measurements after the ambiguity is removed. This is especially 
true for the case they call the implicit ambiguity. 


5.3. Estimate of the CME column density 

We And that the measurement of column density obtained from 
white-light images generally matches the actual column density 
of the MHD simulation to within 10%. The measurements are 
carried out by assuming that the entire column of plasma is lo- 
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Fig. 19. Idealized cartoon showing the location of the real centre of 
mass (red dashed line) and of the folded density centre of mass (orange 
dashed line) for different orientations of the CME with respect to the 
POS (see text). 


cated at one location along the LOS, either on the POS (as usu¬ 
ally done in data analysis), or at the position of the centre of 
mass as inferred from the polarization ratio technique. Our study 
shows that when we assume that all the plasma is located at the 
position inferred from the polarization ratio we get a more pre¬ 
cise measurement of the column density. This is particularly true 
near the flux rope location which is the more interesting region 
for estimating the mass of the CME. In our study the error re¬ 
duces by a factor of 3 when the LOS assumption is used. Hence, 
this can be considered an improved technique for the determina¬ 
tion of column densities from white-light images of CMEs. The 
uncertainties can be minimized by assuming that in each pixel 
of the 2D coronagraphic image the emitting plasma is centred 
on the position along the LOS inferred from the polarization ra¬ 
tio technique. 

5.4. Time evolution 

Einally, when we apply the comparison to a set of snapshots from 
the MHD simulation, we first need to identify the lines of sight 
that cross the CME in order to apply the polarization technique 
to a sensible region of the EOV and to extract information such 
as its trajectory or expansion rate. The polarization technique 
can show several dynamic features, but not all of them belong 
to a coherent and structured ejection like a CME. Provided that 
the relevant lines of sight are identified, the polarization ratio 
technique is able to reproduce the trajectory of the CME and 
its extension with the same spatial approximation that we have 
found when focusing on a single frame. In our study, the offset 
between the polarization technique cloud and the centre of mass 
cloud leads to an error in the CME trajectory of 5° in the longitu¬ 
dinal coordinate. This is a significant discrepancy that could be 
reduced with more complex data analysis that takes into account 
the position on the solar disk where the CME originates. When 
the initial offset is known, it is possible to infer the correct trajec¬ 
tory. The error introduced on the speed of the CME is more dif¬ 
ficult to infer because it is closely connected with the detection 
technique that determines the position and thus the velocity of 
the CME front. On the other hand, the shape of the CME cloud 
and its evolution is captured by the polarization technique. The 
polarization ratio technique enhances a more accurate measure- 
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ment of mass involved in a CME. Especially by using the LOS 
assumption the technique returns a column density accurate to 
within 2% when a CME crosses the field of view. 

5.5. Application of the technique for space weather forecast 

As the polarization ratio technique is able to provide accurate 
measurements of the mass involved in a CME (through the col¬ 
umn density estimate) and the position and velocity of the CME 
(through the centre of mass at different times), we can imagine 
using this technique to detect the early stages of CMEs directed 
towards the Earth and with the support of space forecasting mod¬ 
els to infer the arrival time and the mass. However, as we out¬ 
lined in Sec |5.2[ the polarization technique introduces some er¬ 
rors in the localization of the CME that need to be reduced so 
that the technique gives the best results. The technique is much 
more accurate at localizing the folded centre of mass than the 
real centre of mass. Also, we have highlighted that this error 
increases the more the CME crosses the POS. 

L4 



Fig. 20. Sketch of the ideal configuration for a satellite to continuously 
monitor the initiation of CMEs from the solar corona and to apply the 
polarization ratio technique to compute the 3D structure and the speed 
of the ejection (not to scale). 


As explained in [Bemporad & Pagaii^ ( |20I5| ) and confirmed 
in this work, the technique generally performs best when the 
CME travels at about 20° off of the plane of sky. This peculiar 
position not only minimizes the error due to plasma crossing the 
plane of sky, but also the error introduced by applying the tech¬ 
nique at larger radial distance. We therefore suggest that an ideal 
setup to timely and accurately detect CMEs and measure their 
position, speed, and mass is to place a satellite that observes a 
large field of view in white light (both total and polarized bright¬ 
ness) on an Earth orbit travelling 70° ahead or behind the Earth, 
as in Eig|^ Within this orbit the direction of propagation of 
CMEs directed towards the Earth makes an angle of 20° with 
the POS. This setup would allow us to continuously monitor the 
portion of the corona that generates Earth direct CMEs and to 
promptly alert advanced space weather forecast systems of the 
injection of a CME into interplanetary space with accurate esti- 
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mation of its mass, position, and speed. It should be noted that 
the LA and L5 Lagrangian points are approximately located at a 
60° angle in front of and behind the Earth and thus quite close 
to the ideal position for the 3D reconstruction of CMEs with the 
polarization ratio technique. This is a significant and additional 
advantage to be considered for future space weather monitoring 
missions that will be located in the L4 or L5 Lagrangian points. 
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